Dynamics of orientational ordering in fluid membranes 
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We study the dynamics of orientational phase ordering in fluid membranes. Through numerical 
simulation we find an unusually slow coarsening of topological texture, which is limited by subdiffu- 
sive propagation of membrane curvature. The growth of the orientational correlation length £ obeys 
a power law £ oc t w with w < 1/4 in the late stage. We also discuss defect profiles and correlation 
patterns in terms of long-range interaction mediated by curvature elasticity. 
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Fluid membranes exhibit a variety of shape transfor- 
mations due to their flexibility and sensitivity to changes 
in temperature and osmotic pressure. They are most sim- 
ply modeled by a structureless deformable surface with 
bending energy and a few constraints on its global geom- 
etry jlj . More diverse and prominent deformations occur 
in membranes with internal degrees of freedom, especially 
when they undergo a transition between different ther- 
modynamic phases. An example is provided by a multi- 
component membrane in two-phase coexistence, defor- 
mation of which is controlled by a local coupling between 
composition and curvature. The equilibrium morphology 
as well as the dynamics of phase separation have been ad- 
dressed in a number of theoretical studies j|^J.Q.||,|U^,^,^] , 
which have relevance to the so-called budding and other 
experimental observations. 

Another class of shape transformations is character- 
ized by a coupling to an orientational degree of free- 
dom. It represents anisotropic in-plane shape and/or 
configuration of the constituent molecules. The possibil- 
ity of quasi-long-range orientational order has been pro- 
posed Jl0| , pdp^ ] and the resulting morphologies (such as 
tubules, sponges, and egg-cartons) and defect energetics 
have been explored |11 i^, |l^,|i^,|l5l,[i"6 Jl7 LA special case 
of this is the nematic ordern]l||l|Jj4||§ , the presence 
of which is suggested in the so-called "gemini" sur- 
factants (each of which is a pair of molecules covalently 
linked at their polar heads) pq |. It is recently shown 
that the nematic order, in combination with chirality, 
accounts for the formation of helical ribbons by the ma- 
terial |§. 

In this paper, we address the nematic phase order- 
ing from a dynamical perspective, which has been lack- 
ing so far. Via numerical simulation, we show that 
the coarsening of topological texture is greatly deceler- 
ated by the coupling to membrane deformation. The 
growth of texture size is asymptotically characterized by 
a power law with an exponent less than 1 /4, instead of the 
1/2 for ordinary nematic liquid crystals |u],|2(],|2l|] . We 
provide the interpretation of this novel non-equilibrium 
effect. To prepare for this, we also study how the 
curvature-orientation coupling controls defects and cor- 
relation structure in mechanical equilibrium. Our results 
follow from a Ginzburg-Landau model in the simplest 
setting of a flat topology under a weak coupling condi- 



tion. It would be straightforward to include more com- 
plicated though realistic aspects, such as the stiffness 
anisotrop y |]To|,^6| , non-flat topologies El ] , and interlayer 
coupling [ pL4|Jl 5| | . 

We consider an initially flat membrane without any 
overhang, and describe its profile in the Monge gauge as 
z = h(r) = h(x,y). The orientational order parameter 
Qij(r) = Qij{x, y) is a symmetric traceless tensor related 
to the scalar order parameter S and nematic director n, 
as Qij = S(riinj — 5ij/2)(i,j — x, y). The model free en- 
ergy consists of three parts, which are the homogeneous, 
Frank, and curvature-elastic contributions. Under the 
assumption |V/i| < 1, we retain the lowest order terms 
with respect to V/i, as 
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(summation over repeated indices is implied). Here, A, 
C, and M are the Landau-de Gennes coefficients |2^|, k 
is the bending rigidity, and a is the coupling constant. 
Finally, EL is the unit tensor for an non-chiral nematic 



membrane 12 and the totally antisymmetric tensor for a 
chiral nematic [[l6| . From now on we concentrate on non- 
chiral membranes, but discussions are parallel for chiral 
ones, as we shall see. Neglecting long-range hydrody- 
namic interactions mediated by the solvent, we write the 
kinetic equations in the form 



dQi 



( SF 



dt v \6Qij 

where (s) denotes the symmetric traceless part, and 



(5) 



|J = = -T h (kV 2 V 2 /i + adidjQij) . (6) 



To analyze the initial growth of order parameter upon 
a quench into the nematic phase, we linearize Eqs.(^|^) 
with the aid of the unitary transformation 
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C = a,Tctan(q y /q x ). 
The three eigenmodes are expressed as 
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and eo(q) = Q-(q), where c(g) = -nq 2 /a + Tq(A + 
Mq 2 )/(T h aq 2 ). Their growth rates 7„ = d(lne a )/dt 
(a = +, — , 0) read 
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and 70(9) = — rQ(A + Mq 2 ). The spinodal point is lo- 
cated at A = a 2 jK HQ as seen from the behavior of 
7 a (g) at small q. 

The nonlinear dynamics of phase ordering is studied 
by numerically integrating Eqs.(|^,^) on a square lattice. 
We choose A = -1, C = 20, M = 1, K = 20, and 
a = 1 as the standard parameter set, with the mesh size 
Air = 1. The Landau coefficients give the equilibrium 
scalar order parameter S eq ~ y / 2\A\/C ~ 0.3 and the 
defect core size l CO re ~ \J M/\A\ = 1. The large ratio 
k/M = 20 is not unrealistic if we consider the large stiff- 
ness k ~ 20fcsT usually found in biological membranes. 
We impose periodic boundary conditions on a 512 x 512 
lattice. For the initial condition, random numbers uni- 
formly distributed in [—0.1, 0.1] are assigned to each com- 
ponent of Qij . The kinetic equations are integrated using 
the Euler scheme with Tq = 0.1, Th = 0.1 and the step 
increment At = 0.025. Note that the flat membrane ap- 
proximation breaks down as the texture size exceeds the 
equilibrium curvature radius ~ K/(aS eq ). All the data 
shown below are taken from the time region in which 
(|V/i|) < 0.15 (we denote a spatial average by (• • •)). 

Shown in Fig.l are snapshots of the Schlieren texture 
Q 2 y (r) and the mean curvature V 2 /i(r). In the late 
stage, the coarsening proceeds via pair annihilation of 
topological defects with s = +1/2 and —1/2, where s 
is the winding number for the apolar vector n (or the 
disclination strength j23|). It is seen that a —1/2 de- 
fect accompanies three lobes of positive mean curvature, 
while a +1/2 defect has one. Where the director points 
to the defect core, V 2 /i is negative (positive) for a +1/2 
(-1/2) defect (resp.). To understand this, we note that, 
the membrane shape in the latest stage of coarsening is 
locally well equilibrated with respect to the order param- 
eter configuration. Then, for analysis of static correlation 
properties, we can assume the mechanical equilibrium 
condition SF/Sh — 0, or 



V 2 (V 2 /i) = didjQi 



(10) 



This can be solved in the Fourier space as h(q) = 
(a/ Kq 2 )Q+{q), substitution of which into (^) gives the 
effective curvature elastic energy as 



^ \Q + {q)t 



(27T) 



(11) 



It becomes dominant over Frank elasticity when the tex- 
ture size is larger than the characteristic length, 



(12) 



Within this distance from a defect, the optimum configu- 
ration Qij(r) is little affected by the elastic coupling and 
approximately minimizes Frank elastic energy. By fur- 
ther assuming S = const., n — (cos 9, sin 9), and the po- 
lar coordinate system (r, </>) , the approximation gives |23|] 



s{4> — </>o), 4>o = const. 



(13) 



for a defect of index s located at the origin. The Poisson 
equation (fLOh for mean curvature has the special solution, 



^ 2 h = 
ab 
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This qualitatively explains the mean curvature profile ob- 
tained by simulation. It also shows that a +1/2 defect 
carries much larger bending energy than its anti-defect. 
On the other hand, the anisotropic part of the curvature 



Hij = didjh - (V 2 h/2)6i 



(15) 



is proportional to Qij in a uniform equilibrated system, 
and the principal curvature axis is parallel to the direc- 
tor. This holds also in the coarsening system in the late 
stage and except in the vicinity of defects, as we con- 
firmed numerically. 

The long-range elastic interaction also modifies ori- 
entational correlation at lengthscales larger than A. In 
the real space, and under the approximation S — const., 
Eq.(pll) is rewritten as 



a 2 S 2 
" ,rv ~ Hi-/; 



drdr 



, cos 2 [9 re i(r,r') + 9 re i(r',r) 



where re i(r, r') is the angle between n(r) and r' — r [ p7[ . 
It suggests that orientational correlation is enhanced in 
directions parallel and perpendicular to the local direc- 
tor, while it is suppressed in oblique directions. This 
can be quantified by the relative orientation correlation 
function, which we define as 

G rel (r - r') = (Q«(r)Qy (r + U(n(r)) ■ (r' - r))) 

where U(n(r)) is a matrix of rotation that maps n(r) 
onto e x . Simulation result in Fig. 2 shows that the ex- 
pected tendency is prominent in the region \r — r'\ A, 
where correlation in oblique directions is negative. 

The coarsening kinetics is monitored through the ori- 
entational correlation length £q, which is defined by the 
correlation function Gq(\t — r'\) = (Qij(r)Qij(r')) via 
Gq(^q)/Gq(0) = 1/2. Similarly, the curvature correla- 
tion length £jj is defined as the half-value decay length 
of G H {\r - r'\) = {H t] {r)H l0 {r')). In Fig.3, the ori- 
entational correlation length is plotted as a function of 
time. In the absence of coupling, £q is well fitted by 
a power law £q oc t w ° with wo = 0.43 ± 0.02, in good 
agreement with previous results PUM . Remarkably, the 



elastic coupling makes the growth much slower. The in- 
stantaneous growth exponent 

w Q (t) = dQji^Q)/d(hxt) (16) 

decreases from an initial value ~ wq, and then converges 
to a much smaller value as t — > oo. Deviation from 
the initial value starts earlier for a stronger coupling. In 
the intermediate stage, WQ(t) shows an undershoot and 
even turns negative for most of the parameters we stud- 
ied. The asymptotic exponent Woo is estimated to be 
0.12 ± 0.02 for a = 2 and 0.13 ± 0.01 for a = 3 @. 
On the other hand, the curvature coarsening exponent 
WH^t) = d(ln£ii)/d(lnt) first decreases monotonically 
and then converges to the same value Wqc within er- 
ror bars. In the bottom of the figure we plot the ratio 
£q/£h, which decays to unity after an initial increase. In 
the same figure, we indicate the time at which the ratio 
i\fcurv\) I (/f) reaches 2. It shows that the two lengths 
merge as the curvature-orientational coupling becomes 
dominant over Frank elasticity. 

How does the coupling decelerate coarsening? First we 
note that the effective elastic interaction (|ll]) does not 
provide an explanation. In the mechanical equilibrium, 
the free energy density f curv is reduced to its homoge- 
neous minimum except in a region of size ~ A around 
each defect, in which Frank and curvature elastic energies 
are balanced. Thus the ^Q-dependence of f cur v is that 
of the defect density oc 1/£q, which is also the scaling 
of Frank elastic energy (except for a logarithmic correc- 
tion). Therefore, the effective elastic interaction cannot 
affect the coarsening exponent. This was confirmed by an 
additional simulation that directly implements the long- 
range interaction using Fourier transformation. It shows 
that (i) coarsening curves £q(£) and are both well 

fitted by the exponent 0.43 for null coupling; (ii) dissi- 
pation rates c?(/f) /dt and d {fcurv) /dt are of the same 
order and their ratio is roughly constant in the late stage. 

Thus we can attribute the origin of slow coarsening to 
the dynamics of shape change, Eq.@. It is subdiffusive 
and the source term (oc didjQij) is small except near de- 
fects. According to a simple dimensional counting, the 
mean defect separation scales like £q and the timescale of 
curvature relaxation is proportional to £q. If Woo > 1/4, 
the coarsening of curvature field cannot keep up with that 
of order parameter, and hence / c „ increases. This is 
energetically disfavored once the coupling becomes dom- 
inant over Frank elasticity. Therefore, the asymptotic 
exponent cannot exceed 1/4. The reason why is even 
smaller than 1/4 is not very simple. A possible explana- 
tion is the logarithmic correction to the ^Q-dependence 
of Frank elastic energy. In the XY model and ordinary 
nematic fluids in 2D, it causes a logarithmic correction to 
the power law p9[ and the apparent exponent is smaller 
than 1/2 in finite time simulations |gj2l|,||,||,|6| . A 
similar effect is expected for the present case, because 
the coarsening is still driven by Frank elasticity. 

The non-monotonicity of £g(t) in the intermediate 
stage corresponds to the fact that Frank elastic energy is 
temporally sacrificed to reduce curvature elastic energy, 
which has initially grown on the smaller scale Be- 
cause £q grows faster than £h in the initial stage, it is at 
least necessary that wq < wh holds temporally for the 
two lengths to finally merge. 



Some remarks are in order. First, hydrodynamic 
flow of the solvent effectively modifies the mobility as 
Fh(q) = 1/(4^0') (rj: viscosity) Q. However, the cur- 
vature dynamics is still subdiffusive, and the coarsening 
exponent cannot exceed 1/3. We have performed a sim- 
ulation incorporating the hydrodynamic interaction, by 
solving Eq.(6) with Th(q) above in the Fourier space and 
transforming back each time step to solve Eq.(||) in the 
real space. We found the coarsening exponent £q to be 
0.17 ± 0.03 (for a — 2 and 77 = 1.25), which confirms the 
weakness of hydrodynamic accerelation. The details of 
this simulation will be presented elsewhere. 

Second, the mechanism of slow coarsening is unique 
to membranes with a continuous degree of freedom and 
is irrelevant to phase separation dynamics, which is de- 
scribed by a scalar order parameter . This has the 
following reason. In phase separation, significant changes 
in the order parameter occur only at domain boundaries. 
However, domain boundaries are also the sources of cur- 
vature propagation, as we can see from the corresponding 
dynamic equation for h. Therefore, the membrane shape 
can immediately follow the composition change. In con- 
trast, the orientational order parameter changes signifi- 
cantly even far from moving defects. In order to adapt 
to this, shape deformations spread from the defects and 
over the whole texture, which is slow. 

Finally we consider the effect of chirality. With the 
antisymmetric coupling between curvature and Qij, a 
straightforward calculation gives the elastic interaction 
as 




The chirality reverses the correlation anisotropy so that 
the director correlation is enhanced in oblique directions. 
However, the coarsening law is unaltered since the last 
term in the membrane's dynamic equation @ is still pro- 
portional to VVQ and sources of curvature propagation 
are localized around defects. 
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FIG. 1. Snapshots of Q^ y {r) on a 128 2 portion of the lat- 
tice, at t = 400, 1600, and 6400. The last panel shows the 
mean curvature V 2 /i(r) with the director field in short lines. 




FIG. 2. Relative orientation correlation function G re i(r) 
at t — 6400 (averaged over 10 independent runs) . 
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FIG. 3. Top: orientational correlation length £q versus 
time. The lines corresponding to wq — 0.13 and 0.43 are 
drawn as guides to the eye. Bottom: the ratio between orien- 
tational and curvature correlation lengths. Arrows indicate 
the time at which {\f C urv\} / (/f) = 2. 



